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Abstract 

We consider production of dilepton pairs through coherent electro-magnetic 
radiation during nuclear collisions. We show that the number of pairs produced 
through bremsstrahlung is about two orders of magnitude smaller than the yield 
measured by the CERES collaboration. Therefore, coherent bremsstrahlung can 
be ruled out as an explanation for the observed enhancement of low mass dileptons 
in CERES and HELIOS. 



Production of low mass dileptons in ultra-relativistic heavy ion collisions is consid- 
ered a useful probe of possible in-medium changes of hadrons due to the onset of chiral 
restoration (see |l| for a recent review). Dilepton are penetrating probes, i.e. once 
produced they do not re-interact with the hadronic environment and thus, provide in- 
formation about the early stages of the collision where high temperatures and densities 
are reached. Since vector mesons such as p, u and $ have direct decay channels into 
dileptons, possible in-medium changes of the vector mesons masses can be observed in 
the dilepton invariant mass spectrum. 

There has been a major surge of interest in dilepton mass spectrum due to a recent 
experimental observation of an enhanced dilepton yield at invariant masses of about 
400MeV in S + Au and Pb + Au collisions in SPS compared to proton induced collisions 
as reported by CERES collaboration M. There is also a similar enhancement reported 
by HELIOS collaboration at more forward rapidities 0. Although the decay of the final 
state hadrons explain the data well in proton induced collisions, they can not explain 
the current SPS data for S + Au and Pb + Au in the mass region around ~ 400 MeV . 

A major source of this enhancement is simply pion annihilation which is not present 
in the proton-nucleus system. Including pion annihilation and contributions from other 
hadronic reactions the calculations reach the lower end of the sum of statistical and 
systematic errors of the CERES data (see [Q for a compilation of different calculations). 
An additional enhancement can be achieved || |] if one assumes that the mass of the 
p is lowered according to the conjectures of |7]]. Another possible enhancement arises 
from an in-medium modified pion dispersion relation. While this effect is small if one 
considers the modification of the pion dispersion in a pion gas it is considerably larger 
if, in addition, one takes the effect of baryons into account ||. Consistency of the latter 
scenario with the observed pion to baryon ratio is presently debated. 

Another source of dilepton production is simply the bremsstrahlung due to the de- 
celeration of the incoming nuclei during the collision. Motivated by a recent result of 



Mishustin et.al. [|T0|], we investigate whether this source could account for part of the 
observed enhancement of dilepton numbers as reported by CERES collaboration. In 
reference Jl0| only the production of uj mesons due to the deceleration and their subse- 
quent decay into dileptons has been taken into account. This corresponds to considering 
the iso-scalar part of the electro-magnetic current. However, one also has to take into 
account the iso-vector part, i.e. production of p-mesons and their subsequent decay into 
dileptons. These two amplitudes, iso-scalar and iso-vector, interfere destructively and as 
a result the dilepton production cross section from the coherent deceleration scales like 



the square of the charge and not the square of the baryon number, as assumed in [10 



We should also note that the point like limit of this process addressed here has been 
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considered in 







As already mentioned, one expects that the coherent radiation will be enhanced by a 
factor of Z\%2 in nuclear as compared to proton induced collisions. Here Z\ and Z2 are 
the atomic numbers of the colliding nuclei. The enhancement factor follows from 
assuming coherent radiation off of charged nuclei. This assumption is an approximation 
which should be valid when one considers photon virtualities (invariant dilepton masses) 
which are much smaller than inverse size of individual nucleons. Therefore, as one goes 
to higher and higher invariant dilepton masses, this approximation will cease to be valid. 
Also, the ratio of photon virtuality to center of mass energy of the collision should be 
small so that there is also longitudinal coherence. Even though incoherent radiation will 
become as important and eventually dominate over coherent radiation as one goes to 
higher and higher masses, the coherent radiation will always be there as a background 
and so therefore, it should be understood. Here we try to provide an upper limit on this 
background. 



1 Radiation from Decelerating Nuclei 

Let us consider a typical ultra-relativistic nuclear collision where nuclei A and B move 
towards each other with very high but constant velocities. In order to simplify the 
extremely complicated process of nuclear collisions, we will assume that the main effect 
of the collision is deceleration of each nucleus while they are passing through each other. 
Since electric charges in the nuclei are decelerated, they emit photons. All time like 
photons with virtuality q 2 > which subsequently decay into dileptons are emitted 
during this time. After a passing time t ~ r a+Rb ^ move on w jth reduced but 

constant velocity. We will ignore all subtleties associated with expansion of nuclei in the 
transverse direction during the collision and, for simplicity, assume a Gaussian form for 
the charge distribution. 

We can relate the number of dileptons produced in this process to the Fourier trans- 
form of the correlator of the electro-magnetic currents of the colliding nuclei. It is given 
by0 



where p^ 1 is the virtual photon momentum, a = 1/137 is the electro-magnetic coupling 
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constant and W„ v is the Fourier transform of the product of the electro-magnetic currents: 



WV(p) = / d A xd A ye- tp{x - y) J^x)Jl(y). 



(2) 



Our task is now simple, we just need to write down an electro-magnetic current 
corresponding to an extended charge density with (proper) time dependent velocity. It 
is @: 



dTv^{r)5[v v {x u - z u {t))} 



1 + a v {x v - z v {t)) 



(3) 



with 



and 



z (r) = z (T) + / dr'j(r') 

JT 



z 3 (t) = z 3 (T) + / dr'j(r')P(r') 

JT 



We choose the initial time T = — oo for convenience. f[(x — z) 2 } is a properly normalized 
but otherwise arbitrary charge profile at this point. 

Here, r is the proper time and a M = -p is the corresponding acceleration (decelera- 
tion) with 



v^{t) = 7(r) 



/ 1 \ 




V P{t) J 



(4) 



being the four velocity and 7(7-) = [1 — /5 2 (r)] -1 / 2 . It is easy to verify that this current 
is conserved. The acceleration term a v looks peculiar and arises out of a consistent 
application of the concept of simultaneity in special relativity to an extended (but still 
rigid) charged object and is essential for current conservation in the case of a spatially 



extended charge distribution (see [13] for a nice illustration of this). It will drop out 
when one takes the point charge limit of this expression. 

Expression @ is the current for a charged nucleus moving from left (z = —00) to 
right [z = +00) with velocity (3. In order to get the current of a charged nucleus moving 
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from right to left, we simply take (3 — > —f3 in the current of the right moving nucleus. 
The total current to be used in @ is the sum of the currents of the right moving and 
left moving nuclei. The last step is to determine the velocity (3. We will assume that 
both nuclei have a constant initial velocity $ until they collide at r = 7$. During the 
collision, from Tj to Tf, velocity changes in a non-trivial way. After r = t/, both nuclei 
have again a constant but reduced velocity (3f. 

To proceed further, we need to take a specific form for the nuclear charge distribution 
f[(x — z) 2 ]. For simplicity we will use a Gaussian profile 



f[(x - zf}= p exp 



(t - ^p(r)) 2 - (g - z(r)f 
2a 2 



(5) 



Here po and a are related to atomic number Z and radius of nucleus R by 

1 



po = Z(27riT/3)"3, a' 



-it* 



A more realistic profile would be a Woods-Saxon shape. However, since we want to 
provide an upper limit of the bremsstrahlung contribution, Gaussian distribution is a 
reasonable approximation. For a given charge, it is more narrow and will lead to a larger 
dilepton yield at finite momentum than the corresponding Woods-Saxon profile. 

In order to model the deceleration phase, one could take the velocity during the 
collision to be a linearly changing function of time but it is perhaps more realistic to 
take the rapidity, rather than the velocity, to be a linearly changing function of time. 
We therefore take the rapidity during the collision to be 



V\J) 



Vi + 

At 



(6) 



where (3{r) = tanh[y(r)], Ay = yj — yi, At = 77 — Tj and the initial and final rapidities 
and times are to be determined by experimental considerations. 

Fourier transforming the current to momentum space, dividing the proper time inter- 
val into three different regions and performing the proper time integration in the initial 
and final stages where the velocity is constant, we get: 



J 3 i? (p) = Zexp 



exp 



ii?HW-p 3 ) 2 



i[p° -/3ip 3 - it] 
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- Pf- 



exp 



1 JPtJC/V - p 3 ) 2 + i%\p\lfPf ~ lib) ~ P\lf ~ 7,)] 



i[p° - {3fp 3 + ie] 



! f T i{t)3(t) 



.R 2 Ay 



x exp 



-^ 2 7 2 (t)[/?(t)p° -p 3 ] 2 + ^[P°(l(r)(3(r) - 7i &) -p 3 ( 7 (r) - 7 ,)] 



(7) 



where we have set Tj = without loss of generality. The first two terms (proportional to 
fa and (3f) correspond to the contribution of the initial and final stages of the collision 
respectively while the last term is the contribution of the time interval when the nuclei 
are passing through each other and has to be evaluated numerically. Also, in the last 
term corresponding to the contribution of the deceleration region, extended structure of 
the source is responsible for the second term (proportional to area of the nucleus, R 2 ) in 
the bracket : 



(8) 



and will be absent in collision of point charges. 

Using current conservation PqJq{p) = p 3 J 3 {p) and also J M = + , we then have 



d 4 N, 



i+i- 



a 



1 



dMdyd 2 p t 6vr 3 M 



1 _d 

2 

Po 



Hp)A(p). 



(9) 



As a consistency check, it is easy to show that one does not get either physical (light 
like) or time like photons without acceleration. 

Below, we plot our results for the number of dileptons produced for two different 
values of collision times with and without the CERES acceptance cuts for Ay = 2.4 14 



From here, it is clear that bremsstrahlung is irrelevant for the observed enhancement of 
dilepton spectrum even with our Gaussian charge distribution which would clearly over 
estimate the number of produced dileptons and that a different mechanism is needed. 
Variation of our parameters like rapidity change Ay and collision time Ar does not 
change our results appreciably even in the extreme case of Ay = y$ i.e. full stoppage of 
nuclei after collision. 

We also plot the transverse momentum spectrum of the produced dileptons, inte- 
grated over rapidity and dilepton invariant mass range [200 — 600 MeV] and compare it 
with CERES data. Both graphs have the peculiarity that as one increases the time it 
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takes the nuclei to pass through each other from 1 fm to 4 fm, the number of dileptons 
increases in contrast to what one expects intuitively (at least for high invariant masses). 
This is again due to the extended structure of the source. Unlike the case of a point 
charge, we have an inherent scale in the problem, namely R, the nuclear radius. 
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Figure 1: Number of dileptons produced for lfm and 4fm collision times with and 
without CERES cuts. 

Physically, there are two competing effects which are responsible for this initial in- 
crease in the number of emitted dileptons as one increases the collision time. The first 
effect is just what one expects; as charges decelerate over a longer time interval, they 
emit fewer dileptons. The second effect is due to the finite time required to build up 
coherence among different parts of the extended source. In other words, if the nucleus 
decelerates too quickly, there is no time for the different points in the extended source 
to communicate and coherent emission takes place only from a fracture of the source 
which has had time to react. As we increase the collision (deceleration) time, we in- 
crease the fraction of the nucleus which can be considered a coherent source and as a 
result, we have more dileptons emitted. After some time, the whole nucleus is emitting 
dileptons coherently after which the first effect takes over and the number of dileptons 
starts decreasing as we increase the collision time any further. 
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Figure 2: P t spectrum of dileptons produced for 1 fm and 4 fm with and without CERES 
cuts integrated over rapidity and mass in the mass range between 200 and 600 MeV. 

It is interesting to note that dilepton mass spectrum due to deceleration of a point 
charge (not shown here) shows a (modulated) periodicity which depends on the time it 
takes the nuclei to pass through each other and in principle could be used to determine 
this time. However, this structure is totally wiped out by the Gaussian charge distribu- 
tion and it is unlikely that it could be experimentally useful in determining the collision 
time for a realistic charge distribution. 

In summary, we consider production of dilepton pairs due to coherent bremsstrahlung 
in ultra-relativistic heavy ion collisions. We provide an estimate of the upper limit for 
this contribution using a Gaussian charge distribution. We find coherent bremsstrahlung 
to be a negligible source for dileptons. 

This work was supported by the Director, Office of Energy Research, Office of High 
Energy and Nuclear Physics Division of the Department of Energy, under contract No. 
DE-AC03-76SF00098 and DE-FG02-87ER40328. 
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